Two-Mode Squeezed States and Entangled States of Two Mechanical Resonators 
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We study a device consisting of a dc-SQUID with two sections of its loop acting as two mechanical 
resonators. An analog of the parametric down-conversion process in quantum optics can be realized 
with this device. We show that a two-mode squeezed state can be generated for two overdamped 
mechanical resonators, where the damping constants of the two mechanical resonators are larger 
than the coupling strengths between the dc-SQUID and the two mechanical resonators. Thus we 
show that entangled states of these two mechanical resonators can be generated. 

PACS numbers: 03.67.Mn, 85.25.Dq 



I. INTRODUCTION 

Motivated by their relevance to quantum information, 
coherent quantum behavior of macroscopic solid-state de- 
vices are of great interest. Quantized energy levels, coher- 
ent time evolution, superposition and entangled states — 
have all been observed in various solid-state devices, such 
as quantum dots and superconducting quantum inter- 
ference devices (SQUIDs). Nanomechanical resonators 
(NAMRs)i*2i2i4 with frequencies as high as Giga Hertzs 
can now be fabricated£&£i£. At milli-Kelvin tempera- 
tures, such mechanical resonators are expected to ex- 
hibit coherent quantum behavior. In order to detect and 
control mechanical resonators, some transducer methods 
must be used. These include optical methods, magneto- 
motive techniques, and couplings to single electron tran- 
sistors)i2j2ii. A novel design of mechanical qubits based 
on buckling nanobars was recently studied in Ref. 9. 
Also, buckled modes analogous to buckled-bars have been 
proposed for magnetic nanostructuresM, and mechanical 
bars bent by electric fields have been considered, e.g., 
Ref. [111. Moreover, these quantum mechanical nanobars 
can exhibit behavior similar to superconducting quantum 
circuits^. 

The quantization of NAMRs has been studied 
by coupling NAMRs to a superconducting charge 
qubits— : 14 i 15 ' 16 i 17 . By controlling the charge qubit, a 
NAMR can be prepared into different quantum states. 
Also, charge qubits can also be used to measure the quan- 
tum states of the NAMRs. Quantum nondemolition mea- 
surements of a NAMR were studied with an rf-SQUID 
acting as a transducer between the NAMR and an LC 
resonator—. It was shown that a strong coupling cavity 
QED regime can be realized for a NAMR and a super- 
conducting flux qubits or a NAMR with magnetic tip 
coupled to an electron spin 2 ^. 

For their usages of sensitive displacement detection be- 
yond the standard quantum limit, single-mode squeezed 
states of a nanomechanical resonator were theoretically 
studied. It was shown that squeezed states of the 



nanomechanical resonator can be generated by either pe- 
riodically flipping a superconducting charge qubit cou- 
pled to ili^i or by measuring the superconducting charge 
qubit coupled to it 22 . The squeezing of the nanoresonator 
state can also be produced by periodically measuring its 
position by a single-electron transistor—. 

Two-mode squeezed states, when these two modes are 
from two spatially-separated macroscopic objects, are 
macroscopic entangled states. The generation of these 
entangled states of macroscopic objects are of funda- 
mental interest. Several protocols have been proposed 
to entangle two tiny mirrors with the assistance of pho- 
ton a 24 : 25 . Here we study a device consisting of a dc- 
SQUID with two opposite sections of the SQUID loop 
suspended from the substrate. The suspended parts, 
shaped as doubly-clamped beams, can be approximated 
as NAMRs. The magnetic flux threading the loop of the 
dc-SQUID is modulated by the displacements of both 
NAMRs. Then the dynamics of the dc-SQUID is modi- 
fied by the NAMRs. We study how the potential energy 
of the dc-SQUID is modified by the displacement of the 
NAMRs. We show that the nonlinear coupling between 
a dc-SQUID and the NAMRs, where the dc-SQUID is 
approximated as a quantum harmonic oscillator, offers a 
flexible method for the detection and control of NAMRs. 
Specifically, we discuss two-mode squeezed states of these 
two NAMRs through an analog of the two-mode paramet- 
ric down- conversion process in quantum optics. We show 
here that two-mode squeezed states of the two NAMRs 
can be obtained even when the couplings between the 
dc-SQUID and the two NAMRs are weaker than their 
damping rates. In contrast to this, for the proposal stud- 
ied in Ref. [26| the entanglement between the NAMRs de- 
creases rapidly with time. Other forms of squeezing of 
mechanical oscillators (phonons) have been studied about 
a decade ago in Refs. |23HIH3l 

This paper is organized as follows. At the beginning 
of Sec. |TT] our device is described. Then, after we study 
the potential energy of the dc-SQUID, the Hamiltonian 
of the device is presented. The interaction Hamiltonian 
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FIG. 1: (Color online) Schematic diagram of the top view 
of the device considered here. It consists of a rectangular- 
shaped dc-SQUID and two mechanical resonators shown in 
blue on the left and right sides. These two opposite segments 
of the dc-SQUID are freely suspended and are treated as two 
nanomechanical resonators (NAMRs). The dotted lines indi- 
cate the equilibrium positions of the left and right NAMRs. 
Each "x"in the top segment of the loop represents a Joseph- 
son junction. Xl (Xr) is the displacement of the center of 
the left (right) resonator. The two NAMRs could be located 
sufficiently "far apart" by using a SQUID with an appropriate 
aspect ratio. This could provide EPR-type correlations on the 
two NAMRs located sufficiently "far apart" for wide-enough 
SQUIDs. 



between a dc-SQUID and two NAMRs is complicated 
and has many terms. However, if the frequency of the 
dc-SQUID is properly chosen by the bias current of the 
dc-SQUID, then only a few terms dominate the dynam- 
ics of the coupled system, which is illustrated by writing 
the interaction Hamiltonian in the interaction picture. 
Then, in Sec. Mil we study a special case where, under 
an appropriate choice of the parameters, the interaction 
Hamiltonian is simplified to study the two-mode para- 
metric down-conversion process in the device. Squeezed 
states of the two NAMRs are studied by the Heisenberg- 
Langevin method. Conclusions are given in Sec. HVl 



II. COUPLING A DC-SQUID WITH TWO 
NANOMECHANICAL RESONATORS 



The device we studied is schematically illustrated in 
Fig. Q] It consists of a dc-current-biased SQUID with 
rectangular shape and with two mechanical resonators. 
The left and right sides of the SQUID are suspended 
from the substrate and form the two mechanical oscil- 
lators, our NAMRs. We assume here that these two 
doubly-clamped beams vibrate in their fundamental flex- 
ural modes and in the plane of the SQUID loop. We 
use the following notations II (Ir) for the current in 
the left (right) Josephson junction, and tfL (<Pr.) for the 
phase drop in the left (right) Josephson junction. The 
two Josephson junctions are assumed to be identical and 
have the same critical current I c . Thus the bias current 



h of the dc-SQUID has the form 

h = I c (sin <p L + sin ip R ) 



(1) 



We assume that the inductance of the dc-SQUID loop 
is negligibly small, and thus the magnetic energy of the 
circulating current in the dc-SQUID loop is neglected. 
Thus the voltage drop over the two junctions is zero. 
Therefore, <PR. — tpL = <Pt, where ip t is the phase related 
to total magnetic flux <!> t threading the dc-SQUID loop 



(fit = 2ir 



(2) 



Here $o = h/2e is flux quantum. Introducing two new 
variables 



(<Pr + <Pl) , 

(<PR - <PL) , 



(3a) 
(3b) 



and taking into account that (p_ 
in Eq. (QJ can be written as 

lb = 2/ c sin(^cos 



ipt/2, the bias current 



(4) 



It is here assumed that Xl (Xr) is the amplitude for the 
fundamental flexural mode of the left (right) beam. Let 
Bl (Br) be the magnetic field normal to the plane of the 
SQUID loop near the left (right) mechanical beam and 
<!>{, the external applied magnetic flux threading perpen- 
dicularly the dc-SQUID loop when X L = X R = 0. It 
is assumed that Bl (Br) is constant in the oscillating 
region of the left (right) beam. Then, the total magnetic 
flux threading the dc-SQUID loop is given by 



(5) 



where <$>x is the additional magnetic flux when the two 
NAMRs are displaced from their equilibrium positions: 



$ x = B L X L l + B R X R l 



(6) 



Here, I is the effective length of the left and right beam. I 
is defined as I = Sl /Xl , where Sl is the area between the 
equilibrium position of the NAMR and its bent configu- 
ration. Namely the area Sl spans the region between the 
blue dashed line and the blue bent line in Fig. Q] Equa- 
tion © indicates that the variables of the two NAMRs 
enter in the dynamics of the dc-SQUID by influencing 
the flux threading the dc-SQUID. The influence of the 
two NAMRs on the dynamics of the SQUID can also be 
revealed quantitatively in the potential energy of the dc- 
SQUID, since this is also a function of the displacements 
of the two NAMRs. Thus, we first study the potential 
energy of the SQUID and afterwards the entire Hamilto- 
nian of the coupled system. 
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FIG. 2: (Color online) The potential energy U(<p, $x) (scaled 
by Ej) of the dc-SQUID as a function of the phase variable ip 
and the magnetic flux <!>x originating from the displacements 
of the two NAMRs. The red color represents higher potential 
energy U , and the blue color represents lower potential energy. 
Both ip and <5>x are shown in units of 7r/$o. In (a)-(c), the 
bias magnetic flux threading the loop of the dc-SQUID is 
set at 2n$o, (2n + J)$o, ano - + h)^o, respectively; and 



the bias currents are all set at lb 
and It — 0.5J C . 



0.1/c 



In (d), $i, = 2n$o 



A. Potential Energy of the Vibrating dc-SQUID 

The potential energy of the dc-SQUID is 



Ufa* 



x 



-2EjCOS 7T 



x 



* 



cos <p - ^-Ejtp , (7) 



where Ej = HI C / (2e) is the Josephson energy of the junc- 
tio n 31 i 32 . To have an idea of the situation under which 
the dc-SQUID can be described by a quasi-particle in a 
quadratic potential, in Fig. [2] (a)-(d) we plot this poten- 
tial energy |J7J| for various values of the bias magnetic 
flux *b and the bias current Since *x/*o < 1 in 
the case of experiments using a GHz NAMR, here we 
focus on the limit *x/*o < 0.1 in Fig. [2l A particle 
in a quadratic potential can be described by a harmonic 
oscillator when its kinetic energy is much smaller than 
the barrier of the potential. We notice that, for the 
bias magnetic flux <&£, = 2n$0i with n being an inte- 
ger, the phase variable ip falls in a potential well when 
the NAMRs oscillate around their equilibrium points. It 
is possible to approximate the dynamics of tp as a har- 



monic oscillator. The charging energy of the dc-SQUID, 
E c = (2e) 2 / (2Cj), is assumed here to be much smaller 
than the modified Josephson energy (cosqo)Ej of the dc- 
SQUID. Here, q = sin" 1 (I b /2I C ) is a value of ip which 
corresponds to one of the minima of the potential energy 
U(ip,*x), when the two NAMRs are at their equilib- 
rium positions. Cj is the capacitance of the left and 
right Josephson junctions. 

We then expand the potential U {ip, $x) near one of its 
minimum points ((p,*x) — (<70)0). When *x/*o "C 1, 
and = 2n$o, the first cosine in Eq. depends weakly 
on <E>x- Up to second order in $x, we can also expand 
After shifting the origin of ip to qo, and omitting the 
constant terms, the potential energy in Eq. becomes 



4> 



Ej(cosq ) tp 2 -Ej(l- cos g ) ( w^- 



E, 



(singo)^ + ^(cosq )p 2 



x 



<t>0 



(8) 



Therefore, if the first term in the above potential is much 
larger than the other two terms, the dynamics of ip will 
still be well described by a harmonic oscillator. 

The higher order terms, such as 7r 4 $ x / (24<&o), in the 
expansion of cos (ir*x/*o) are negligibly small for the 
situation considered in our paper, when *x/*o ~ 10~ 3 . 
Theoretically, it is possible to increase the ratio 7r<I>x / 
by using a stronger magnetic field Bl and Br and/or 
using soft NAMRs with greater zero-point fluctuations. 
However, in practice, the magnetic field is limited to (up- 
most) in tens Tesla and the zero-point fluctuations of the 
NAMRs is less than 10~ 12 m for the most of the exper- 
iments. Therefore, the periodic nature of the Josephson 
Hamiltonian has no chance to play a role here. Indeed 
the situation considered here is very similar to the op- 
tical parametric-down-conversion system, except for the 
coefficients of polynomial expansions of the interaction 
Hamiltonian. 

Now we consider how well a dc-SQUID is approximated 
by a harmonic oscillator. Since the barrier of the poten- 
tial U (tp,*x) has a finite height, the dynamics of p is 
not an ideal harmonic oscillator. However, if the energy 
of the quasi-particle is small enough, the dynamics of ip 
can be approximately described by a harmonic oscillator. 
The maximum number N max of energy levels that can 
be confined in the potential U (ip, 0) is A^ max = AU/fl, 
where the height of the potential is 

AU = 2Ej [2 cos q + sin q (2q - %)] , (9) 

and Q is the frequency of the harmonic oscillator. 



B. Hamiltonian of the Coupled System 

Near the minimum potential U(qo,0), the free Hamil- 
tonian of the dc-SQUID can be written as a harmonic 
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oscillator Hamiltonian 



where 



E c ^ 2 + E'jv 2 , 



with 



E'j = Ej cos q , 



(10) 



(11) 



where the constant term has been omitted when we de- 
rived the Hamiltonian in Eq. (|10p . It is convenient to 
introduce the annihilation and creation operators a and 



,.t. 



E' 



1/4 



Eli V + 



E r 



K 

E' 



1/4 



1/4 



(12a) 
(12b) 



Then the free Hamiltonian in Eq. ([TOj) of the dc-SQUID 
can be rewritten in the form 



H s = Ml a 1 a 
with the angular frequency 



ft — \ E r E', 



J 



(13) 



(14) 



When the energy of the dc-SQUID is not very large 
{{a) a) < -/V max ), the dynamics of ip is well described 
by a harmonic oscillator under a suitable bias magnetic 
flux threading the loop of the dc-SQUID. It is conve- 
nient to also introduce annihilation and creation op- 
erators for the fundamental flcxural modes of the two 
NAMRs (i = L, R) 



bi = \l — Xj + i 



2h 



t = 

2h 



1 



■Pi 



X; - i 



1 



IhvriiLOi 



P 



(15a) 
(15b) 



Here, Xi and Pi are the coordinate and momentum oper- 
ators of the ith NAMR; nij and u>i are the effective mass 
and angular frequency of the ith NAMR. The effective 
angular frequency u>i is not the one of the fundamental 
flexural mode, which is modified by the second term in 
the potential Eq. ijHJ). Then the free Hamiltonian of the 
two NAMRs can be written in the form 



H 



NAMR 



Mo L b\b L + Mo R b\b R , 



(16) 



where the constant terms have been omitted. Thus, 
in terms of creation and annihilation operators, from 
Eq. I|8|). the interaction Hamiltonian between these two 
NAMRs and the dc-SQUID are given by 



V 



9L 



9R 



ci a 



(a + a f ) +c 2 (a + a^Y 



(17) 



9l 



9R 



ci 



C-2 



$o V 2m L ujL 



2m R uj R 



O / Ej cos q 

-(tango) 

n 



1/4 



(18a) 
(18b) 

(18c) 
(18d) 



The interaction Hamiltonian (fiTf is central to this 
work. Notice that it contains both linear and nonlin- 
ear terms. Generally, it is very difficult to evaluate the 
behavior of this coupled system. However, since the fre- 
quency O of the dc-SQUID can be set by the bias cur- 
rent h, we can reduce the interaction Hamiltonian V 
in Eq. (|17|) to a simplified form by invoking the rotat- 
ing wave approximation. We now rewrite the interaction 
Hamiltonian V in Eq. (|17p in the interaction picture, with 
the free Hamiltonian 



Hq = Mia) c 



hio L b\b L 



f)ijj R b R b R 



(19) 



Then the terms of the interaction Hamiltonian V can be 
classified by the ways that the frequencies u>l, Wr, and 
fi can be combined. In Table Q] we list half of the cou- 
pling terms and the combinations of their frequencies. 
The other half are their corresponding Hermitian conju- 
gate terms, which have the same frequencies but with a 
negative sign. 

In Table [H it can be seen that, for large detuning, one 
needs to mainly consider the zero-frequency terms in the 
first row of the table. Then this interaction Hamiltonian 

V enables a quantum nondemolition measurement of dis- 
crete Fock states of a NAMR, as discussed in Ref. 0. 
When the frequency of the dc-SQUID is set at some spe- 
cial value, one can mainly consider the resonant terms. 
For example, if the frequency of the dc-SQUID and those 
of the two NAMRs are properly set so that Vi ^ ojl = u> R 
and also f2 ^ 2lul = 2u> R , then only the zero- frequency 
terms and resonant terms in the interaction Hamiltonian 

V are kept under the rotating wave approximation. The 
reduced interaction Hamiltonian V r consists of the terms 



V r = ci9l9r b R bb a) a + H.c. 



(20) 



which in fact offers us a mechanism for coupling two 
NAMRs. Thus, our proposed device offers a flexible 
(literately) model for the control and measurement of 
NAMRs. 



III. TWO-MODE SQUEEZED STATES OF TWO 
NANOMECHANICAL RESONATORS 

In this section we focus on the two-mode squeezed 
states of the two NAMRs. It is possible to produce 
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Frequencies 


Interaction terms 








T ~ ^ I ^ I ~ T 

2 C2(g L b L b L + gRb R b R )a i a, 


1 


2lo l 


C2gib 2 L a t a, 


2 


2lor 


C2g R b R a'a, 


3 


LOL + LOR 


C2gLgRbLbRa ] a, 


4 


LO L - LOR 


c 2 gLgRb R bLa < a, 


5 


20. 


C2{gib\b L +g R b R b R )a 2 , 


6 





ci(glb L b L + g 2 R b R b R )a, 


7 


2lo l +20. 


C2g L b L a , 


8 


2lo l - 20 


c 2 g 2 L b 2 L a t2 , 


9 


2lo r + 20 


C2g 2 Rb 2 Ra 2 , 


10 


2lor - 20 


c 2 g 2 R b 2 Ra i2 , 


11 


LO L +LOR+ 20 


c 2 gLgRbLbRa 2 , 


12 


LOL + LOR — 20 


c 2 gLgRbLbRa) 2 , 


13 


LOL — LOR, + 20 


c 2 gLgRbLb R a 2 , 


14 


LOL — LOR — 20 


c 2 gLgR.bLb R a} 2 , 


15 


2lo l + 


cig 2 L b 2 L a, 


16 


2lo l - 


Cl glbla ] , 


17 


2lo r + 


Cig 2 Rb 2 R a, 




ZLOR — i I 


^,2 u2 „t 

cigRb R a' , 


19 


LOL + LOR +0 


cigLgRbLbRd, 


20 


LOL + LOR — 


agLgRbLbRo) , 


21 


LOL — LOR. + 


cigLgRbLb R a, 


22 


LOl — LOR — 


cigLgRbLb R a\ 



TABLE I: Terms in the interaction Hamiltonian Eq. (117[) and 
their frequencies, in the interaction picture. 

entangled states of the two NAMRs by considering the 
analog of the parametric down-conversion in quantum 
optics. The zero-frequency terms in Table U commute 
with the free Hamiltonian (fT3")) of the dc-SQUID and the 
free Hamiltonian (fT(I|) of the two NAMRs. Let us assume 
that the proposed circuit works at low temperature. If 
the two NAMRs are initially in the vacuum state or in 
very low-energy states, then we have 5lr <SC fi, with 

S L R = c 2 (g 2 R {b\ i b R )+gl{blb L )) . (21) 

Then we can rewrite the free Hamiltonians of the dc- 
SQUID Eq. ([HI) and the two NAMRs Eq. HU) as 

H' a = Q'a^a + io L b[b L + uj R b R b R , (22) 

where Q' = £1 — 5lr- By properly setting the bias current 
lb one can let tt' ~ lo l — lo r — 0. Then, in the interac- 
tion picture, after adopting the rotating wave approxi- 
mation, we simplify the interaction Hamiltonian between 
two NAMRs and dc-SQUID as 

V = i 1 (a)b L b R + ab\b R ^ , (23) 

where 

V = -ci9l9r ■ (24) 

Driven by this interaction Hamiltonian V' , two-mode 
squeezed states of the two NAMRs can be produced in 



the device similarly to a light beam interacting inside a 
nonlinear medium in quantum optics, because both of 
them follow the same Hamiltonian (|23p . 

We now consider that the mode of the dc-SQUID is 
in a coherent state \a), where |cv| 3> 1. Then we can 
treat the mode of the dc-SQUID as a classical field and 
replace the operator a in the Hamiltonian V in Eq. |23|) 
by a complex number |a| exp (— i<p). Then, in the inter- 
action picture defined by the Hamiltonians (f22j) and (|23|) , 
the dynamics of the coupled system is described by the 
following Hamiltonian 

Vi = e* \a\ V b L b R + e~** \a\ r, b[ b R . (25) 

The motions of b L and b R are 

b L (t) = cosh(7)6i -ie^sinh(7)6jj , (26a) 
b R (i) = cosh (7) b R - ie~^ sinh (7) b\ , (26b) 

in the interaction picture of the Hamiltonians (|22|) and 
(j25| . with 

1 =\a\r ] t. (27) 

The generation of two-mode squeezed states of these 
two NAMRs can be shown by their collective coordinate 
and momentum operators 

X T [t) = X L (t)+X R (t), (28a) 
P T {t) = P L {t)+P R {t), (28b) 

where, Xi (t) and Pi (t), i = L, R, are defined by Eq. (fT5|) 
by substituting hi and b\ with bi(t) and b\(t) in Eq. (|26|) . 
The uncertainty relation for the collective coordinate and 
momentum operators Xt (t) and Pt (t) is 

A [X T (t)} A [P T (t)} = h |cosh 2 7 + e 2 ^ sinh 2 7 | . (29) 

In Eq. (f2"9"]) we have assumed that the zero-point fluctu- 
ation of positions of the left NAMR 

S L = V / /i/(2m LWL ), (30) 
and that of the right NAMR 

S R = y/h/(2m R uj R ) (31) 
are the same. Here 

5 X = V26 L = V2S R (32) 

is defined as the zero-point fluctuation of the collec- 
tive coordinates Xt of the two NAMRs. And (p — 
= V%Cr is defined as the zero-point fluctuation of 
the collective momentums Pt of the two NAMRs. Here, 
C 2 = V(2<5 2 ), i = L,i?. 

If we choose cj) — — tt/2, then the variance of the col- 
lective coordinates A>r(i) becomes 

A [X T (t)} = 5 X exp( 7 ) . (33) 



Notice that 7 < because 7 = — c\gLgR.\ct\ t. Therefore, 
perfect two- mode squeezed states, i.e., pure entangled 
states, of the two NAMRs are generated. 

The variance of Xx{t) (the entanglement) was ob- 
tained above by assuming that both the left and right 
NAMRs be in their ground states. It can be checked that 
if both the left and right NAMRs are initially in coherent 
states or thermal states, then the Hamiltonian (|25[) will 
not produce entangled states of them. However, if only 
one of the NAMRs is initially prepared into a number 
state, then entangled states of these two NAMRs can be 
generated by the Hamiltonian ([2"5)l . For example, when 
the left and the right NAMRs are initially prepared in the 
number states |0) and |1), respectively, then the Bell-type 
entangled state ai|01) + a2|10) can be generated. Here, 
cti and 02 are complex numbers. When one of the two 
NAMRs is initially in a coherent state and the other one 
is in the vacuum state, then the so-called "single-photon- 
added coherent states"— can be generated by the Hamil- 
tonian ([25]) . 

Let us now consider the more realistic case where both 
the dc-SQUID and the two NAMRs are coupled to their 
environments. The quality factors of the two NAMRs 
with GHz frequency are smaller than that of the dc- 
SQUID 6 ' 34 . The quality factor of a GHz NAMR is of 
the order of 10 3 , while that of a superconducting circuit 
can be as large as 10 6 . Therefore, below we consider the 
noise from the environment acting on the two NAMRs. 
To include damping effects, due to the noise from the 
environments, on the dynamics of the two NAMRs, we 
use the Heisenberg-Langevin equation method 3 -^. Then, 
for the motions of the operators of the NAMRs, we have 
the following set of equations: 



dt 
dt 



Fl (t) , 
F R (t) , 



As in the ideal case in Eq. ([55)) we also let 
Here, 

£ = \a\v 



(34a) 
(34b) 
-tt/2. 

(35) 



is the effective coupling strength between the dc-SQUID 
and two NAMRs. Also, Kl and kr represent the damp- 
ing rates of the left and right NAMRs, respectively; and 
the associated noise operators are Fl (t) and Fr (t). We 
evaluate the properties of the states of the two NAMRs 
by the variance of the collective coordinates Xt ■ We find 
that the damping of the two NAMRs help producing two- 
mode squeezed states of the two NAMRs, regardless of 
the initial states. The variance of the collective coordi- 
nates Xt is calculated as 



[A(X T )] 2 = 



s 2 s 2 

-±- (2k r A c + «_) + (2k l Af - k-) 
-8^£A 5 , 

K + 



(36) 
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FIG. 3: (Color online) The squared variance of the collective 
coordinates Xt of the two NAMRs as the function of damping 
rates kl and kr of the left and right NAMRs, both kl and 
kr are normalized by the effective coupling constant £. Here, 
Sx is the zero fluctuation of the collective coordinates of the 
two NAMRs. 



under the Markov approximation and in the overdamped 
case: £ < kl/2 and £ < Kr/2. In the Appendix A we 
outline the main ideas of the derivation. Here, 



K± = K L ± Kr , A ? = 



klkr - 4£ 2 



(37) 



When the zero-point fluctuations of the left and right 
NAMRs are equal, we have 



[A(A T )] 2 = f Acfl-— ^— 

4 V K L + Kr 



(38) 



Fig. [21 shows, the variance of the collective coordinates 
Xt versus the damping rates kl and kr. It is clear 
from Fig. [3] that appreciable squeezing can be gener- 
ated even when the dampings of the two NAMRs are 
severe (ten times the coupling constant £). This indicates 
that the squeezing is robust against damping. The max- 
imum squeezing is obtained when both damping rates 
(for the left and right NAMRs) approach the coupling 
strength between them and the dc-SQUID. Since the cou- 
pling strength £ is proportion to \a\, one can increase the 
squeezing rate by gradually increasing the power of the 
microwave applied to the dc-SQUID. As the damping 
rates of the two NAMRs increase, the squeezing effect 
decreases steadily. 

To consider the experimental feasibility of our pro- 
posal, we choose the following parameters for the two 
NAMRs and the dc-SQUID 



JUL 



I = 



kl 
E c 
Ej 



wiR = 1CT 18 
1.5 GHz, 
1.2 GHz, 

10 /ilTl, 

Br = 1 T, 
Kr = 2 MHz 
0.061 GHz, 
120 GHz. 



Kg, 



(39a) 
(39b) 
(39c) 
(39d) 
(39e) 
(39f) 
(39g) 
(39h) 
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It was already demonstrated in experiments^^ that a 
10 /um long doubly-clamped beam can oscillate with a 
frequency of several Giga Hertzs. The effective mass 
of this antenna-shaped beam is much smaller than its 
weight. And its effective mass can be further modified 
when the beams are under strains and stresses. The 
numbers used here for the dc-SQUID are also consis- 
tent with the experimental numbers shown in Ref. l37l . 
Then the additional magnetic flux from the two NAMRs, 
$x/^o — 5 x 10~ 3 <C 1, would satisfy our assumption in 
section [TTJ The maximum number -/V max of energy levels 
confined in the current biased potential energy is calcu- 
lated as iV max ~ 150. Therefore, the harmonic oscillator 
approximation and the classical field approximation for 
the dc-SQUID are both possible. The time needed to 
obtain the two-mode squeezed state is determined by the 
effective coupling constant £. Assuming the same damp- 
ing rates for the two NAMRs, the maximum squeezing 
A (Xt) = (1/2) Sx can be obtained. Therefore, it should 
be possible to realize our proposal of generating two- 
mode squeezed states of the two NAMRs with current 
experimental conditions. 

To experimentally detect the generated two-mode 
squeezed state of the two NAMRs, a (in principle) rela- 
tively direct method would be checking the variance of 
the collective coordinate Xt of the two NAMRs. Since 
the left and right NAMRs are symmetric in the interac- 
tion Hamiltonian (|17p , they can be treated as one virtual 
NAMR. To detect two-mode squeezed states of the two 
NAMRs, the detection methods should be able to ap- 
proach standard quantum limit of the NAMRs. With tra- 
ditional displacement detection methods ^ 39 ' 40 ! 41 , such 
as optical interferences, magnetic-motive method and 
coupled single electron transistor, the best record of de- 
tection precision was about 4.3 standard quantum lim- 
its^. There are also other proposals for displacement 
detection by coupling the mechanical oscillator to some 
two level syste m 13 i 19 i 21 . These methods in principle can 
detect quantum states of the NAMRs. After the entan- 
gled state is generated, one can switch the dc-SQUID 
to the phase qubit regime 4 ^, and utilize the nonlinear 
coupling between the virtual NAMR and the dc-SQUID. 
Then the SQUID can be used to measure the variance of 
Xt, as discussed in Ref. [U 



netic flux through the dc-SQUID, thereby these modes 
cannot be coupled to the dc-SQUID, even through their 
frequencies match the resonant condition. It is similar to 
the case of the experiments of magnetomotive detection 
of flexural oscillation of NAMRsi, where torsional and 
strain-stress oscillations have been neglected. 

As mentioned in Section. IIIH to generate two-mode 
squeezed states of NAMRs, the NAMRs should start 
in their ground states or number states. The NAMRs 
should be cooled such that the thermal excitation energy 
is less than those corresponding to the NAMRs' frequen- 
cies. For a one-GHz NAMR, this means that the temper- 
ature should be below 50 mK, which is still within the 
capability of dilution refrigerators. In principle it is pos- 
sible to prepare the NAMRs used in our proposal in their 
ground states. Moreover, recently there have been many 
efforts in reducing the temperature of mechanical res- 
onators by active coolin g 43 ' 44 i 45 ' 46 i 47 . Also a temperature 
as low as 5 mK was already demonstrated for a mechan- 
ical resonator—. Besides, there are also theoretical pro- 
posals for the production of number states of NMARs 48 . 
Therefore, though currently the ground states and (or) 
number states of NAMRs might be difficult to prepare 
experimentally, we expect these to be realized more ear- 
lier in the near future. For example, if the temperature 
of a GHz NAMR reaches 10 mk, the thermal occupation 
number will be ~ 10~ 3 , which is nearly a ground state. 

In conclusion, we have proposed a device to couple a 
dc-SQUID to two NAMRs, which can be used to create 
an effective coupling between these two NAMRs, and also 
to measure and control the two NAMRs. We have shown 
that two-mode squeezed states can be generated in a ro- 
bust fashion by this device, in analogy to the two-mode 
parametric down-conversion process in quantum optics. 
This two-mode down-conversion process offers us a pro- 
tocol of producing entanglement in two mechanical res- 
onators in a solid state device, while previous proposals, 
see, e.g., Refs. |24||25||49| . were based on entanglement- 
swapping by the assistance of photons. Our proposal 
might be promising for the experimental test of the exis- 
tence of entangled states of macroscopic objects. 



IV. DISCUSSIONS AND CONCLUSIONS 

In our proposal, we only consider the fundamental vi- 
bration modes of the NAMRs. Generally, there are also 
vibration modes with higher frequencies, torsional and 
strain-stress oscillations in the NAMRs 1 -. The vibration 
modes with higher frequencies will be excited only when 
they happen to resonate with the dc-SQUID, which can 
be easily avoided by optimizing the parameters of the 
NAMRs. As for torsional and strain-stress oscillations, 
they are hardly coupled to the dc-SQUID. These modes 
of oscillations of the NAMRs will not change the mag- 
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APPENDIX A: HEISENBERG-LANGEVIN 
EQUATION FOR TWO NANOMECHANICAL 
RESONATORS 

We consider the Markov approximation (Fl^ r ) = 
F[ r \ = 0. Using Eqs. ([55kV([Mb). a solution of the 
expectation values (&£ (t)} and (b R (t)} can be given as 



(in) = 



b L (0) cosh (ft) - b R (0) sinh (ft) 



(Ala) 



b R (0) cosh (ft) - b[ (0) sinh (ft) 



(Alb) 



It is seen that below the thresholds f < kl/2 and f < 
kr/2 we have 



(6l> = (b R ) = 



(A2) 



The variance of the collective coordinates Xt can be eval- 
uated by the expectation values of the bilinear operators 
of the two NAMRs. These are the expectation values of 
the quadratic operators of the left NAMR 



Li = (bl) , L 3 = (#) 



L 2 = (b{b L + b L bl 



(A3a) 
(A3b) 



the expectation values of the quadratic operators of the 
right NAMR 



Ri 
R2 



b% 



{bl) , Rs 



b R b R + b R b R ) , 



(A4a) 
(A4b) 



and the expectation values of the quadratic operators of 
both NAMRs 



Ci = (b L b R + b R b L 
C 2 = 



b L b R 



b R b L ) = C) 



t 

4 i 
ft 

3 • 



(A5a) 
(A5b) 



From Eq. (f3~4")) it is found that these expectation val- 
ues satisfy a closed set of equations of motion^. To de- 
termine the values involving the expectation values of 
the products of the noise operators and the operators of 
the NAMRs, we rewrite Eq. (|3"4"a V (f3~£b ) and their corre- 
sponding Hermitian ones in the matrix form 



where B 



B = -MB + T , 



b L (t),b[ (t),b R (t),b R (t) 



(A6) 



and T 



represents the transpose operation. Here 



r «l n 
2 u 



M = 





f 

A formal solution of Eq. 



f 
f 




Kt 
2 

2 

4f 



Bit) 



B(0) 



6|) is given by 



(A7) 



(A8) 



Multiplying the above equation by J 7 ' (i) from the right 
side, we obtain 



B (t) F (t) 



t u\ - 



-Mt 



B{0)F* (t) 



- M (*-*)F(1?)dlfjA (t) 



(A9) 



Since the operators of the NAMRs at the initial time 
t = are statistically independent of the noise opera- 
tors, we have (0(0) F (t)) = 0. Using the fact that the 
corresponding elements of the matrix of the left part of 
Eq. (|A9[) and those of the matrix of the right part of 
Eq. (|A9[) are equal, and combining the Markov approxi- 
mation, we obtain 



M*) *!(*)) = % 

bR (t)F R (t) 



2 

K R 

2 



(AlOa) 
(AlOb) 



All other products of the operators of the two NAMRs 
and the noise operators are zero. Therefore, in the in- 
teraction picture, finally when the expectation values of 
these bilinear operators do not change with time, the so- 
lution of the above set of equations reads 



L\ — L 3 — Ri — R 3 — C 2 — C3 







(All) 



and 



L 2 


k_ 

K+ 


% K R a 

— A c , 


(Al2a) 


R 2 


K- 
K+ 


2«l a 


(A12b) 




= c 4 = 




(A12c) 



with f < kl/2 and f < k r /2. Also, in the interaction 
picture, we have (&l (f)) = {b R (t)) — after a sufficiently 
long time. Therefore, the variance of the collective coor- 
dinate Xt becomes 

[A (A T )] 2 = S 2 L L 2 + 6 2 R R 2 + 5 L 5 R (d + C 4 ) (A13) 



Fl (i) , F\ (t) , F R (t) , F R (t) are vectors, and [...] T This provides the main result of section IIIH 
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